The heating history of Vesta and the onset of differentiation 
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Abstract 

In this work we study the link between the evo- 
lution of the internal structure of Vesta and ther- 
mal heating due to ^^Al and ^°Fe and long-lived 
radionuclides, taking into account the chemical dif- 
ferentiation of the body and the affinity of ^^ Al with 
silicates. 

We explored several thermal and structural sce- 
narios differing in the available strength of energy 
due to the radiogenic heating and in the post- 
sintering macroporosity. By comparing them with 
the data supplied by the HEDs and the Dawn 
NASA mission, we use our results to constrain the 
accretion and differentiation time as well as the 
physical properties of the core. 

Differentiation takes place in all scenarios in 
which Vesta completes its accretion in less than 1.4 
Ma after the injection of ^^Al into the Solar Neb- 
ula. In all those scenarios where Vesta completes 
its formation in less than 1 Ma from the injection 
of ^^Al, the degree of silicate melting reaches 100 
vol.% throughout the whole asteroid. If Vesta com- 
pleted its formation between 1 and 1.4 Ma after 
^^Al injection, the degree of siHcate melting exceeds 
50 vol.% over the whole asteroid but reaches 100 
vol.% only in the hottest, outermost part of the 
mantle in all scenarios where the porosity is lower 
than 5 vol.%. If the formation of Vesta occurred 
later than 1.5 Ma after the injection of ^^Al, the 
degree of silicate melting is always lower than 50 
vol.% and is limited only to a small region of the 
asteroid. 

The radiation at the surface dominates the evo- 



lution of the crust which ranges in thickness from 8 
to about 30 km after 5 Ma: a layer about 3-20 km 
thick is composed of primitive unmelted chondritic 
material while a layer of about 5-10 km is eucritic. 
Keywords Asteroid Vesta; Asteroids, differentia- 
tion; Thermal modeling. 



1 Introduction 

Is Vesta a terrestrial world? And if so, of which 
kind? Vesta is the most geologically diverse of the 
large asteroids, showing a surface geology as differ- 
entiated as the one of the Moon and Mars as re- 
vealed by resolved images taken by the the Hubble 



Space Telesco pe ( |Gaffey||1997[ [Binzel et al.||1997 
Li et al.||2010l and recently confirmed by the spec 



troscopic measurements of the Dawn mission (De 
Sanctis et al.||2012 ). Surface spectroscopy of Vesta 
indicates regions that are basaltic (see, e.g., Mc- 



Cord et al 119701 |Gaffey|1997| [Binzel et al.|1997[ |Li 



et al.|2010[]De Sanctis et al.|2012[ ), therefore due to 
ancient lava flows and effusive phenomena. This is 
surprising evidence that Vesta once had a molten 
interior, like the Earth does. After the crust for- 
mation, Vesta underwent an intense and extended 



resurfacing due to coUisional evolution (Turrini et 



al. 2011 Coradini et al. 2011), as can also be in- 



ferred by the presence of a large crater on the south 
pole ( [Schenket a"L]|2012[ ). 

From a spectroscopic point of view, most of the 
terrains of Vesta have a typical pyroxene signature 



(McCord et al. 1970 De Sanctis et al. 2012), as 
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revealed by the presence of the 0.9 and 1.9 ^m 
absorption bands. The presence and the shapes 



of these bands in the spectra of Vesta, moreover, 
match those in the spectra of the Howardite, Eu- 
crite, Diogenite meteorites (HEDs in the fohowing, 
see, e.g.,|McCord et aL|1970l|Gafrey|1997[|De Sanc- 



tis et ah 2012), a match that led to the considera- 



tion of Vesta as the source of the HEDs cohected 



on the Earth (see McSween et al. 2011 for a more 



detailed discussion) . This would make Vesta one of 
few intact bodies in the Solar System that differ- 
entiated and for which rock samples are available. 
Therefore Vesta can be used as a case study to in- 
vestigate the primordial stages of the evolution of 
the terrestrial planets. 

The differentiation of asteroids requires a strong 
heat source, capable of producing high tempera- 
tures in bodies with high surface-to- volume ratios. 
The heat source must have also operated quite early 
in the history of the Solar System so that the differ- 
entiation could have occurred within the first few 
miUions years of the Solar System (see, e.g., Scott 



2007) 



Several heat sources have been proposed: the 
only one that currently seems to be able to meet 
the requirements, is the heat produced by the de- 
cay of short-lived radioactive isotopes. Among the 
possible isotopes the most effective both in terms 
of concentration and of the production of energy 
seems to be ^""Al (|Urey|l955"l). |Bizzarroet al.|(|2005[) 



suggested that the accretion of chondritic parent 
bodies ^-^1.5 Ma after the formation of Calcium- 
Aluminum-rich Inclusions (CAIs in the following) 
is the limit to initiate the melting. 

A wide range of thermal models of planetesimals 
with ^^Al as the heat source exists in the litera- 
ture: we briefly analyze the most significant ones, 
describing first the models concerning Vesta and 
subsequently those regarding generic asteroids. 

The first study on the differentiation of asteroid 
4 Vesta with a numerical code was presented by 
Ghosh & McSween (1998). These authors demon- 



strated that it is possible to sustain partial melt 
on Vesta for ~100 Ma assuming radiogenic heating 
by ^^Al. They examined the influence of the delay 
time of accretion with respect to the formation time 
of the CAIs and they assumed instantaneous accre- 
tion. They showed that melting and core formation 
did not occur if the time of accretion is more than 3 
Ma after CAIs formation, based on the assumption 
that the eucrites were generated with no more than 
25 vol.% partial melting of sihcate ( Stolper|[l977 ). 



They did not use a fixed temperature at the sur- 
face, where the heat loss is governed by the energy 
balance between the radiation at the surface and 
the interior heating due to the decay of short ra- 
dionuclides. They did not examine the evolution of 
the internal structure such as the percolation of the 
iron through the silicatic matrix and consequent 
core formation. 



Gupta & Sahijpal ( 2010 ) analyzed the differenti- 



ation path of Vesta and of other differentiated as- 
teroids of sizes 20-270 km, assuming a linear ac- 
cretion rate and introducing parameterized convec- 
tion in the molten core and magma ocean. They 
included sintering but they used a constant tem- 
perature surface boundary. The authors explored 
two different scenarios dealing with the origin of 
the basaltic achondrites: in the first the origin is 
linked to partial silicate melting, in the second to 
the residual melt left after the crystallization in a 
cooling magma ocean. The core-mantle differenti- 
ation in both the scenarios initiated subsequent to 
40 vol.% silicate melting. 

In literature we also find many papers dedicated 
to the study of accretion and differentiation of 
generic asteroids. 



Merk et al. ( 2002 ) examined the influence of the 



accretion process on the thermal evolution. They 
concluded that the accretion rate must be consid- 
ered as long as the accretion time was not much 
smaller than the half-life of the radionuclide under 
study. 



Hevey & Sanders (2006) also incorporated con- 



vection in their thermal evolution models of plan- 
etesimals assuming that mantle regions where the 
degree of partial melting exceeded 50 vol.%, i.e. 
a magma ocean (Taylor et al. 1993), would 
converting 



be 



They simulate mantle convection by 
assuming that at temperatures greater than 1725 
K the thermal conductivity increases by three or- 
ders of magnitude. As a result, heat is transferred 
rapidly from the converting interior into the over- 
lying rigid, partially molten zone. These authors 
modeled the primordial history of the planetesi- 
mals by analyzing the sintering process when the 
temperature reaches the value of ^700 K. The au- 
thors concluded that parent bodies of differenti- 
ated meteorites accreted within 2 Ma from CAIs 
before most chondritic parent bodies and they sug- 
gested that molten planetesimals may be a source 
for chondrule melt droplets. 



Sahijpal et al. (2007) performed numerical sim- 



ulations of the differentiation of planetesimals un- 
dergoing a linear accretion growth with both ^^Al 
and ^°Fe as the heat sources. They studied in par- 
ticular the dependence of the growth rate of the Fe, 
Ni-FeS core on the onset time of planetesimal ac- 
cretion (relative to CAIs formation), the (constant) 
accretion rate, the final size of the planetesimal, 
and the ^°Fe/^^Fe initial ratio. They did not use 
a radiation boundary condition but they included 
the sintering and linear accretion. The authors de- 
veloped two scenarios: in the first, the core grew 
gradually before silicate began to melt and in the 
second the core segregated once the silicate had be- 
come 40% molten. In their model they showed that 
the rate and time of core and crust formation de- 
pend strongly on the onset time of the planetesimal 
accretion: significant melting occurred when accre- 
tion was completed before 2 Ma. 
Recently, 



authors showed that the differentiation of planetes- 
imals depends principally on the accretion time and 
metal and silicate segregation are almost simulta- 
neous processes and last between 0.4 and 10 Ma. 
As in IGhosh & McSweenl (|1998|), IMerk et al. 



(2002), Moskovitz & Gaidos (12011) and Sramek et 



al. (2012 1 we do not consider sintering and prefer to 



analyze the post-sintering phase, exploring differ- 
ent cases of porosities. This is due to the large un- 
certainties associated to the assumed initial porosi- 
ties and the corresponding thermal conductivities. 
The initial value (0.001 W m~^K~^) assumed by 
Hevey fc Sanders] ( |2006[ ) and [Sahijpal et al.| ( [2007| 



is based on laboratory measurements of the lunar 
regolith in a vacuum ( Fountain fc West]|1970 l and, 
as pointed by Moskovitz & Gaidos (20111, we do 



not consider lunar regolith as an ideal analog for an 
unsintered planetesimal. Similarly we do not model 



how the migration of silicate melt and in particular 
the redistribution of ^^Al from the interior into a 
crustal layer would affect the thermal evolution of 
planetesimal. In their model, core formation would 
require a bulk melting degree of 50 vol.%. They 
only considered the case of instantaneous accretion 
but they did not include convection or sintering 
and they used a fixed temperature at the surface. 
The authors concluded that differentiation would 
be most likely for planetesimals larger than 20 km 
in diameter and that accreted within approximately 
2.7 Ma from CAIs. 



Moskovitz & Gaidos (2011) studied & Sanders (20061 and Moskovitz & Gaidos (2011). 



accretion, as in Ghosh & McSween (1998), 



As pointed by Moskovitz & Gaidos (2011), recent 



Hevey 



dynamical studies treating the turbulent concen- 
tration of small particles in proto-planetary disks 



(Johansen et al. 2007 2009 Cuzzi et al. 2010) 



Sramek et al. (2012) presented a multiphase 



model for differentiation of planetesimals which 
takes into account phase separations by compaction 
driven melt migration. According to these authors, 
for the bodies with radii greater than 500 km, im- 
pacts provide an additional heat source. They 
included accretion, sintering, radiation boundary 
condition, analyzing the evolution of the internal 
structure. The authors concluded that, depending 
on the accretion rate, melting may start either in 
the center, or directly under the surface, or both: as 



showed that planetesimals can grow "nearly instan- 
taneously" in < 100 yr to sizes of 100 km or larger. 
These timescales are much smaller than the half-life 
of the radionuclides here considered so the accretion 
time can be considered negligible. We also neglect 
the heat released by the accretion process since it 
is insignificant relative the short-lived radionuclide 
contribution (in the unrealistic case in which all the 
kinetic energy is converted in heat, the maximum 
temperature increase possible is about 90 K). As in 



Neumann et al. (2012) and Sramek et al. (2012) we 



noted by Merk et al. ( 2002 ) , accretion phase plays a 



use a Darcy's law model for segregation of metal 
and silicates, which is a distinct advance over the 
earlier Vesta-specific models of Ghosh fc McSween] 
(119981 ) and |Gupta fc SahijpalldmOD . The affinity of 
^Alwith silicates causes, in our model, a different 
distribution of the heat in the evolving structure: 
as a consequence temperatures reach the maximum 
value in the mantle enriched in silicates, as we ob- 
serve in the thermal profiles of [Ghosh fc McSween] 



crucial role in the thermal history of planetesimals. (1998) but not in Gupta fc Sahijpal (2010). We 



Neumann et al. (2012) focused on the different!- also assume a radiation boundary condition, unlike 



ation of small planetesimals (< 120 km) for melt 
fractions smaller than 50 vol.%. They combined 
the calculation of conduction, accretion, sintering, 
melting and melt segregation by porous flow. The 



of Gupta fc Sahijpal (20101, which led to the study 



of the evolution of the crust and the depiction of 
a reliable heating history of Vesta. As |Merk et al.| 
( 2002 ) we include Stefan law formulation to incor- 



porate the effect of latent heat into the evolution. 
In the current work we assume as a critical pa- 
rameter the delay time At^, which incorporates 
both the delay in the injection of ^^Al in the So- 
lar Nebula and the temporal interval covered by 
the accretion time of Vesta. The delay time Ai^ is 
important because it determines the content of ra- 
dioactive material at the onset of the thermal evo- 
lution of Vesta. The longer the delay time, the 
weaker is the available source. We also develop sub- 
scenarios characterized by different post-sintering 
porosity in order to improve the treatment of the 
internal structure of the rocky asteroid, in particu- 
lar the formation and the evolution of the core, by 
studying its physical and chemical properties, and 
to depict the primordial history of 4 Vesta (and in 
general of all rocky asteroid partially or completely 
differentiated), constraining the accretion and dif- 
ferentiation times. This work provides a theoreti- 
cal support to data analysis for the Dawn mission, 
which addresses the investigation of the Vesta's in- 
ternal structure as well as the composition of the 
crust and of the underlying mantle in the region of 
giant impact basins. 



2 Geology and Geochemistry 

Vesta formed at the same time and in the same 
Solar System environment as the terrestrial plan- 
ets, but the asteroid never grew to planetary size 
because of the formation of Jupiter at 5.2 AU, 
which caused disruptive resonances in the Main 
Belt. Recent papers investigated the bombardment 
on Vesta due to Jupiter's formation and to its mi- 
gration ([Turrini et "aI1|2011[ |2012[). Their results 



Markowsky et al. (2005) confirmed that metal 



clearly indicated that the formation of the giant 
planet caused an intense early bombardment in the 
orbital region of the Main Belt. In the most ex- 
treme scenario, for Jupiter migrating by 1 AU, the 
basaltic crust of Vesta would not have survived in- 



tact to the Jovian Early Bombardment (Turrini et 
al. 20121. In the other studied cases, Vesta's sur- 



face is saturated by craters as big as 100 km and 



its crust would have significantly eroded (jTurrini et 
aL]|2011[ [Turrini||submitted[ ) 



Evidence for rapid iron-silicate differentiation of 
asteroids, i.e. the formation of the core, comes 
from is^Hf-^^^W concentration variations in iron 



silicate differentiation must have occurred on the 
parent body of magmatic irons no later than ~1 
Ma after CAIs. Early differentiation is also sup- 
ported by the presence of excesses of ^^Mg from 
the decay of extinct ^^Al in angrites (Bizzarro et 



al.||2001[ ) 



al. 2005) and in the eucrites Piplia Kalan (Srini- 
vasan et alT]|1999[ ) and Asuka 881394 ( [Nyquiest et 



Using i82jjf_i82'y^ ^^^g^ £qj. eucrites, Kleine et al. 



( 2004 ) showed that mantle differentiation in Vesta 



(assumed the parent body of eucrites) occurred at 
4563. 2±1. 4 Ma and they suggested that core for- 
mation took place 0.9±0.3 Ma before mantle dif- 
ferentiation. 

Moreover, the compositions of HEDs have been 
used to estimate the asteroid's bulk composition 
and core mass (jMcSween et al. 2011 Zuber et al. 



2011). The evidence for the formation of a metal- 



lic core is provided by the siderophile elements de- 
pletions in HEDs relative to chondrites. |Righer cS] 



Drake (1997) estimated the mass of the core in the 



range 5-25% of the total mass, based on the models 
of siderophile elements abundances in chondrites; 



otherwise, for Ruzicka et al. ( 1997 ) the mass of the 



core is in the range 4-30% of the total mass, based 
on mixtures of HED minerals that produce chon- 
dritic element ratios and using a bulk density of 
~3540 kgm~^, slightly different with res pect to the 
most recent value (3456 kgm~^) given by Russell et 



al. (2012). Ruzicka et al. (1997) estimated a metal 



lie core of Vesta of <130 km in radius. 

The olivine may be a significant component of 
the mantle, based on the occurrence of a limited 
number of olivine-bearing diogenites and a dunite 



having HEDs oxygen isotopic composition (Mc- 
Sween et al.|2011 ). Ruzicka et al. (1997) suggested 
an olivine-rich mantle ~65-220 km thick. 

The crust of Vesta consists of some mixture of 
basalts (eucrites) and a significant fraction of ultra- 



mafic cumulates (diogenite). Warren et al. (2009) 



estimated the mixing ratio of basaltic eucrite to dio- 
genite in howardite regolith breccias to be constant 
at 2:1. The normative mineralogy of that compo- 
sition is ^70 vol.% of pyroxene and ^25 vol.% of 
plagioclase, corresponding to a value for d ensity of 
~3170 kgm~^. This ratio, as suggested by Warren 



et al. (2009), could be an upper limit (and the real 



meteorites (Horan et al. 1989 Kleine et al. 2004) 



crustal proportions may be closer to 1:1) because 
of bias in the impact sampling of deeper lithologies. 



Ruzicka et al. ( 1997[ ) suggested a lower crustal unit 
(~12-43 km thick) composed of pyroxenite from 
which diogenites were derived and an upper crustal 
unit (~23-42 km thick) from which eucrites origi- 
nated. 

While we know that Vesta is a differentiated 
body, the exact mechanism for the formation of its 
core is still a matter of debate. In fact, in literature, 
there are two schools of thought regarding core sep- 
aration from silicatic matrix. The first scenario for 
the formation of the core is that the melt fraction of 
the silicates is required to be larger than about 50 
vol.% ( |Taylor|1992"{ [Taylor et al.|1993[ ), arguing for 
the presence of an early magma ocean in a planetes- 
imal to form a core. This assumption is supported 
by experimental studies that partial melting of me- 



teorites does not show metal migration (Takahashi 



1983 Walker & Agee 1988). The second scenario 



suggests that iron segregation and possibly core for- 
mation can start already for small melt fractions of 



iron ( Larimer]! 1995 Hewins fc Newsoin||1988[ ) even 
before silicate starts to melt. This assumption is 
supported by the observations of Fe, Ni-FeS veins 



in the acapulcoite-lodanite parent body ( McCoy et 



al. 1997) and by recent experiments suggesting an 



interconnected melt network for pressures below 2- 
3 GPa ( |Terasakiet"aI][2008l ). 



3 The Model 

3.1 Initial and Boundary Conditions 
of the Model 

Our model is based on ID heat conduction with 
radiogenic heat source and with black-body radia- 
tion at the surface. We assume Vesta as a spherical 
body of fixed radius equal to 270 km and initially 
composed of a homogenous mixture of two com- 
ponents, the first one (~25 vol.%) generically re- 
ferred to as metals (essentially Fe and Fe-S) and 
the second (~75 vol.%) generically referred to as 
silicates. This composition is similar to those of the 
H and L classes of the ordinary chondrites, which 



contain significant amounts of metals (McSween 



1990) even if the inferred composition for Vesta 



is slightly different, as it appears to be strongly 



depleted in sodium and potassium (Consolmagno 
&: Drake]|1977 ). The post-sintering porosity ranges 
from 1 to 5 vol.%. We opt for these values of poros- 



ity because the lower bulk density estimate derived 
from Dawn suggests 5 to 6 vol.% in the mantle and 
crust, which is consistent with Vesta's intense bom- 



bardment history as revealed by Dawn (jRussell et 
aI][20T2| 



The initial temperature (Tq) of the body (that 
is also the local temperature of the Solar Nebula) 
is fixed to 200 K ( |Lewis|[T974| : a change of To to 
300 K does not affect the results in any significant 
way. During the thermal evolution, the silicatic 
and metallic fractions per unit of volume change as 
a consequence of differentiation. The physical pa- 
rameters density, specific heat and thermal conduc- 
tivity of each volume unit also change accordingly. 
We impose a radiation boundary condition at the 
surface and a Neumann boundary condition (heat 
flux equal to zero) at the center, as expressed by 
the following equations: 



T(r,t = 0)=ro, 



'dT' 

dr 


surj 


'dT' 
dr 



= 0, 



r=0 

ea 

'k 



\ surf 



n 



(1) 



(2) 



(3) 



where Tgurf is the temperature of the surface, e is 
the emissivity and a is the Stefan-Boltzmann con- 
stant (see Tabic fl]). 



3.2 Physical 
Model 



Description of the 



To numerically study the thermal evolution of 
Vesta we take into account the heating due to de- 
cay of ^^Al, ^°Fe and long-lived radionuclides (e.g. 
238u^235u)^ The initial concentrations of the two 
short-lived radioactive elements ^^Al and ^"Fe to- 
gether with their half-lives are reported in Table 

HI 

The equation of heat transfer in a porous 
medium, assuming local thermal equilibrium so 
that Tgoi = Tiiq = T (here sol stands for solid and 
liq for fluid taking averages over an unit of volume), 
becomes (following Nield fc Bejan|[2006| : 



(pc), 



dT 
'Ik 



K„,\/Tj 



H, 



(4) 



where 



{pc)^ = {l-<t'){pc),,i + <jy{pc) 



liq 



K„ 



(1-0) Ksoi + (t>Kug 



(5) 



the eutectic Fe-FeS system, and sihcate melt gen- 
eration is assumed to initiate at 1425 K (see Ta- 
bic II]). The entire latent heat for melting is as- 
sumed to be expended in a temperature "window" 



are the overall heat capacity and the overall ther- 
mal conductivity respectively; H is the overall heat 
production per unit volume of the medium. The 
surface temperature (Tsurf) is controlled by the ra- 
diation boundary condition (see the eq.pl). The 
volumetric radiogenic heating rate, due to^^Al de- 
cay, following Castillo- Rogez et al. (20071, can be 
expressed as: 



between solidus and liquidus (Ghosh & McSween 
[1998 ). This simplification of temperature "win- 
dows" does not change the whole thermal history 
by being the exact latent heat supplied to cause sil- 
icate and metallic melting. The percolation of the 
metals through the silicate matrix, starting at 1213 
K, is governed by the advection equation. If we as- 
sume Y to be the concentration of the metals, the 
equation reads: 



HAi=pCsr['-^Ai\^H*e 



-xt 



(6) 



dY 
~dt 



v-Vr = 0, 



(10) 



where p is the mean density, Csi is the mass frac- 
tion of the silicates, P^AlJo is the initial concentra- 
tion of ^^Al in kg for kg of silicates, A = ln{2)/TAi 
is the decay constant and H* is the specific power 
production (see Castillo- Rogez et "aL]|2007 and Ta- 
bic [l]). An analogous formulation holds for "^"Fe ( |2006[ ), can be expressed: 
and long-lived radionuclides (for decay information 
see Table of |Castillo-Rogez et al.|[2007[ ). 

Once the melting temperature of Fe-FeS (or of 
silicates) is reached, partial or complete melting oc- 
curs depending on the parameter 



in which v is the migration velocity of the molten 
metal, given by Darcy's law. The chemical diffusion 
is assumed negligible. 

The migration velocity of molten metal, following 



Yoshino et al. (2003 2004) and Senshu & Matsui 



Kr 



-g^p 



(11) 



r-T„ 



Ti 



iq 



T, 



sol 



(7) 



following Merk et al. (|2002|) and assuming a lin- 

The 
and 



where Kjj is the permeability of the silicate 
medium, p, — 0.005 Pa-s is the viscosity of molten 
iron, g is the gravitational acceleration, Ap is the 
density contrast between molten metal and solid 
silicate (see Table II]) . The permeability-porosity 
relationship is expressed by: 



ear growth of x with raising temperature (T) . 
values of the temperature for the initial {Tgoi) 
complete (Tuq) melting temperature of metals (or 
of silicates) are reported in Table II] Also follow- 
ing lMerk~eF'ar (2002), the specific heat (that we 



Kd = 



P 



(12) 



assume not to depend on the temperature) is mod- 
ified through the Stefan coefficient. 



where 4> is the porosity, Vg = 10^"^ m is the grain 
size, /3 = 200 is a geometrical constant and n = 2 
is predicted in an isotropic model with regular pore 
network along the edge of grains with tetrakaidec- 



Ste 



Ldx 
cdT 



L 1 

C Tuq - 



(8) 



ahedron form (Yoshino et al 



As noted by Yoshino et al. 



2004) 



( 2003 ) , if permeable 



to take into account in a simple way the latent heat 
during phase transition: 



(l + Ste). 



(9) 



As in Ghosh & McSween ( |1998 ), iron melting is 
initiated at 1213 K, the melting temperature of 



flow is established, segregation velocities (which 
vary in the range 1-100 m/a) are rapid in compari- 
son with the timescale of core formation predicted 
from the ^^^Hf-^^^W isotope system. During the 
percolation the volume fraction of the metallic com- 
ponent (Y) goes down enriching the forming core 
region. When the temperature reaches the value 
of the 50 vol.% melting temperature of the silicate 



(assumed to be 1725 K, Taylor 1992), the separa- 
tion of two melts occurs and the silicate component 
(X) moves upwards to the mantle region dragging 
^^Al with it. At the end of this phase, the core 
becomes pure metallic because the melted metals, 
being more dense than silicatic ones, sinks to the 
center. During the evolution, the concentration of 
Al grows in the mantle underneath the lithosphere 
while the density profile varies due to the differen- 
tiation and the moment of inertia factor decreases 
starting from the initial characteristic value of 0.4 
for a uniform sphere. 

We investigate several evolutive scenarios, vary- 
ing the strength of the radiogenic sources through 
the delay time (At) and the porosity. The com- 
bined solution of the eqs.Q and (10) led us to 
study the evolution of the internal structure, con- 
straining formation time, size and mass of the core 
of Vesta, the size of the "crust" and the tempera- 
ture profile as a function of the distance from the 
center and of the time. 



3.3 Numerical Procedure 

The numerical solution of the system of differen- 
tial equations Q and ( [lO| is obtained using a ID 
finite difference method (Forward-Time Central- 
Space, FTCS) in radial direction. A spatial grid of 
Ar = 300 m is used. To avoid numerical stability 
problems due to the instability of the FTCS scheme 
we use the Lax correction ( Press et al.||2007 ). 

To ensure the stability of our numerical ap- 
proach, we use an adaptive time increment accord- 
ing to the Courant-Friedrichs-Lewy stability condi- 
tions for each of the physical processes (heat dif- 
fusion, metal percolation, radiation boundary con- 
dition) we consider in our work. Following [Toksoz 
& Solomon ( 1973 ), thermal conduction imposes the 



Finally, while the percolation of metals is taking 



following critical time step: 



At 



cond 






(13) 



In analogy with Toksoz & Solomon (1973), we can 



define the following critical time step associated to 
the radiation boundary condition: 



^^rad — 



(pc)„Arr, 



surf 



"( 



surf 



7^4 

^0 



(14) 



place, we need to solve also eq.(lO) and introduce 
a third critical time step: 



Atr, 



Ar 



(15) 



where v is the velocity of the metal percolation. 
The Courant-Friedrichs-Lewy stability condition 
requires that the time step used in our model sat- 
isfies the following criterion: 



At < min (Aicond, ^trad, Ai 



perc) 



(16) 



Therefore, at each temporal iteration of the pro- 
gram we select the actual time step based on the 
minimum critical time step among those we com- 
puted. As the stability condition requires the ac- 
tual time step to be lower than the critical one, 
we chose to use a value equal to 90% of the small- 
est critical time step as a compromise between the 
competing needs for stability and performances. 
So, our time step is defined as: 



At = 0.9 X min [At^ond, Atrad, At 



perc J 



(17) 



4 Results 

The results of our geophysical and thermal model 
depend on all the various physical parameters that 
we reported in Table [l] and described previously. 
In the following discussion we will first focus on 
the dependence of the evolution of Vesta from the 
delay time (Ai^). The delay time is essentially an 
unknown parameter, but it is critical in determin- 
ing the initial overall abundance of short-lived ra- 
dioactive elements, i.e. the intensity of the source 
of energy, the maximum temperature reached dur- 
ing the evolution and then the cooling behavior of 
the object. The chosen values of Ai^ sample a long 
temporal interval and therefore very different inten- 
sities of the radioactive sources. We investigated 
seven scenarios, labeled as N0-N6: as shown in Ta- 
ble |2] we considered values for Ai^ equal to 0, 0.5, 
1, 1.5, 2, 2.5 and 3 times the half-life of ^^Al, cor- 
responding to values in the range 0-2.15 Ma. 

Before proceeding, we want to stress that the de- 
lay time (Aid) does not coincide with the accre- 
tion time of Vesta. The former represents the time 
between the injection of ^^Al in the Solar Nebula 
(which in principle may not coincide with the time 



Quantity 



Value 



Unit 



Reference 



Final primordial radius (i?) 

Density of metal (solid) (p,net,aoi) 

Density of silicate (solid) {psii,soi) 

Density of metal (liquid) {pmet.Uq) 

Density of silicate (liquid) [psii,iiq) 

Specific heat of metal (solid) {c^et,soi) 

Specific heat of metal (liquid) {craet.iiq) 

Specific heat of silicate(solid) [Csu^soi) 

Specific heat of silicate (liquid) [Cf^nug) 

Latent heat of metal (L^^t) 

Latent heat of silicate {Lsu) 

Metal soUdus (TJJf ) 

Metal liquidus (Tj^,'") 

Sihcate soUdus (T^'Ji) 

Sihcate liquidus (T/>') 

50 vol.% of silicate melting temperature (T^q) 

Thermal conductivity of metal (K^et) 

Thermal conductivity of silicate {Ksu) 

Initial metal volume fraction {Y) 

Initial silicate volume fraction {X) 

Post-sintering porosity {(p) 

Temperature of Solar Nebula (Tq) 

Stefan-Boltzmann constant (cr) 

Emissivity (s) 

Half-life of '^^Al (tai) 

Heat production of "^^Al 

Initial isotopic abundance of ^^Al in ordinary chondrites ([^*'A/]o) 

Half-life of '^°Fe (rpe) 

Heat production of ^^Fe 

Initial isotopic abundance of ^^^Fe in ordinary chondrites ([^'^Fe]o) 



270x10-* 

6300 

3200 

6200 

2900 

600 

2000 

720 

720 

270 

400 

1213 

1233 

1425 

1850 

1725 

50 

3 

25% 

75% 

1% - 5% 

200 

5.67x10-** 

1 

0.717 

0.355 

6.20 X 10"'' 

2.62 

0.068 - 0.074 

(22.5-225) X 10" 



Kg 77i"' 
Kg m"' 
Kg m"' 
Kg m^' 
JKg-^K 
JKg-^K 
JKg-^K 
JKg-^K 
KJ Kg- 
K.J Kg- 
K 
K 
K 
K 
K 
W m-i K 
W m-i K 



Ghosh & McSween|1998 
T7 



Neumann et 



INeum ann et 

Neumann et 

Neumann et 

Sahijpai et 

Sahijpai et 

Saiiijpal et 

.Saliijpai et 

^losh & McS 

Kjliosh & McS 

Ghosh & McS 

.Ghosh & McS 



al. 21 
al. 20T7 

al. 'im 
al. ■lun 

al. 20^^ 
al. 20D7 
al. 2UU7 
al. 20T17 _ 
ween 1998 
ween 199!? 
ween 199f? 
ween.l99ff 



I'aylor 



I'aylo 



TWI 



TWI 



Hevey &: Sanders 200 6 



bramel^ et al. 



bramel^ et al. 



im 



lun 



K 




Lewis 


1974 




W m '^ K ^ 








Ma 


CastiUo-Rogez et al.'2009 


WKg-^ 


CastiUo-Kogez et al. 2009 


ppb 


Castillo-Hogez et al. 2009 


Ma 


Kugel et al. 2U09 


WKg-^ 


Castillo-Koge 
Castiilo-Koge 


z et al. 2007 


ppb 


z et al. 2007 



Table 1: Physical parameter values used in this work. 



of the condensation of CAIs) and incorporates the 
latter, which however is always assumed negligible 
respect to Atd- 

Our first analysis of the results obtained for the 
different scenarios was based on the compatibility 
of the simulated evolution of Vesta with the con- 
strains supplied by the HEDs. As we can see in 
Figsjl] [2] and [3J in N0-N2 scenarios, which imply 
delay times less than 1 Ma, the complete melting 
of silicates is achieved across the whole of the as- 
teroid. In N3 scenario, where the delay time is 1 
Ma, complete melting of silicates is achieved in a 
limited region of the mantle of Vesta but in the rest 
of the asteroid the degree of melting is larger than 
50 vol.% with the only exception of the case where 
the porosity is equal to 5 vol.%. Similarly, in N4 
case where Atd = 1.4 Ma, the degree of melting is 
generally lower than 50 vol.%, except possibly in 
the limited region of the mantle. Finally, if Atd > 
1.5 Ma, silicate melting is either not possible or is 
limited only to a small region of the asteroid. 

The N0-N3 scenarios are compatible with the re- 



sults of Greenwood et al. ( 2005 1 which link the for 



mation of eucrite and diogenite to large scale (> 50 
vol.%) melting of the silicates. This would imply 
Vesta formed no later than 1 Ma after injection 
of ^^Al. It is be noted that in those scenarios in 
which the differentiation takes place, the melting 
of the silicatic component begins in the first 1 Ma 
and the differentiation is completed in about 3 Ma. 
The formation of Vesta between 1 and 1.5 Ma after 
the injection of ^^Al in the Solar Nebula is com- 
patible with the formation of HEDs if eucrite and 
diogenites can form from a partial melt ranging to 
25-50 vol.%. Those cases where Vesta appeared in 
the Solar Nebula after more than 1.5 Ma since the 
injection of ^^Al are not compatible with petroge- 
nesis of HEDs. The conditions to start the forma- 
tion of eucrites and diogenites are always obtained 
within 1 to 2 Ma from the formation of Vesta. In 
Figsfll [2] and l3] we set the maximum temperature 
to 4000 K to ease the comparison: we refer the 
reader to Figs|4]^e),[5]je) and[6]je) for the maximum 
temperature reached in the different cases. 

We observe that the maximum temperature 
reached in the N0-N2 far exceeds the liquidus sili- 
cates melting temperature (1850 K): this happens 
because in our model we do not take into account 
other cooling mechanisms (convection and effusive 
phenomena) other than the conduction and the ra- 



diation at the surface. In Figs|4j [5j [6] for three 
different initial post-sintering porosities, we report 
in (a), (b), (c), (d) and (e) temperature profiles at 
different times (0.1, 0.5, 1.5, 3 and 5 Ma), in (f) 
the maximum temperature vs time profile. In (a), 
(b), (c), (d) and (e) the horizontal hues represent: 
the windows for the melting of metals (green and 
red) and silicates (cyan and magenta). Note that 
we report in Figs. l4]-[6] (in (a) - (e)) the profiles 
of N0-N4 scenarios while we neglect N5 and N6 be- 
cause in these scenarios the differentiation of the 
asteroid does not take place. The general trend 
we observe in the different cases is the following. 
The first phase is characterized by a homogeneous 
heating of the asteroid. The second one is an in- 
crease of the temperature in the region in which 
the metal is depleting and migrating to the center 
of the asteroid. The final phase is the formation of 
the metallic core followed by the migration of the 
silicates (with the ^^Al) towards the surface and 
the increase of the temperature in the mantle of 
the asteroid. When the differentiation occurs, the 
formation of a pure metallic core is possible. For a 
delay of about 1.5 Ma, the formation of the metal- 
lic core is possible only if the porosity is lower than 
5 vol.%. 

We point out that, with the assumed chondritic 
composition and in the limits of ID modeling, in 
those scenarios in which the differentiation takes 
place, the formation of a core of about 60 km (with 
a density of about 6200 kgm~^) is possible. Its 
mass represents about the 2 vol.% of the total mass, 
slightly lower than the minimum value (4 vol.%) 
given by [Ruzicka et al" (19971, while the moment 
of inertia {Mol) is 0.33. The core size (60 km) is 
lower (about 1/2) than Dawn's current estimates. 
As we can see in [Russell et al. (2012), assuming a 
core density of 7100 and 7800 kgm"'', a core size of 
107-113 km and a core mass fraction of about 18% 
were derived using a two-layer mass-balance model 
that can reproduce the gravitational moment J2. 
The difference between Dawn's estimate of the core 
size and the current work is due to the 1-D geome- 
try of our model. In such geometry all radial cells 
have the same unity volume, while in 3-D geome- 
try the contribution of external shells is larger. As- 
suming an initial metal fraction of about 25 vol.%, 
the maximum size of a pure metallic core is about 
60 km. Probably, assuming a more iron-rich ini- 
tial composition, the core would be larger and the 



moment of inertia would be different. 

In FigQa) Vesta heats up maintaining an al- 
most uniform temperature due to the initially ho- 
mogeneous distribution of ^^ Al, in all the scenarios. 
Due to the strong heating source, after 0.5 Ma (see 
FigMb)), very high temperatures are reached in NO 
and Nl and this trend continues after 1.5 Ma (see 
FigWc)) and also involves N2 scenario. In N3 and 
N4 silicate melting is still partial. After 3 Ma (see 
Figlll^d)), only in N4 the complete melting of sili- 
cates is not reached and the temperature profile is 
inside the melting temperature window of silicates. 
In FigQe), after 5 Ma, the general trend for all the 
profiles is similar to that after 3 Ma. We report in 
FigMf) the maximum temperature vs time and we 
observe that the differentiation is possible only in 
N0-N3 scenarios, for which the separation of sili- 
cates and metals occurs. The formation of the core 
takes place between 0.45 and 3.06 Ma (in N0-N3). 

We can observe in Figsjs] and [6] (in (a) - (e)) 
that the evolution of the temperature in the first 5 
Ma from CAIs is quite similar to the previous case 
(i.e. porosity of 1 vol.%). In the case of porosity 
of 2 vol.%, after 5 Ma Vesta is in the cooling phase 
in N3 and N4, while if the porosity is equal to 5 
vol.% the asteroid enters the cooling phase in all 
scenarios, except NO and Nl. 

In our simulations we observed that the outer- 
most region of Vesta never reaches the melting tem- 
peratures of either silicates or metals. This region 
would maintain the original composition of the pre- 
differentiated Vesta and we therefore named it the 
chondritic crust. Alongside the chondritic crust, 
we can derive from our simulations the thickness 
of the outer layer of material that is below the 
melting temperature of silicates (i.e. solid, SL in 
the following). Before the differentiation, the solid 
layer instantaneously coincides with the (shrink- 
ing) chondritic crust but, after differentiation oc- 
curs, it can be interpreted as the sum of the chon- 
dritic crust plus the solidifying eucritic and dio- 
genitic layers. The evolution of the SL (chondritic 
crust) is identified, in the plots of Figs|4J [5]and|6] 
by the outermost intersection between the temper- 
ature profile and the line corresponding to the onset 
of silicate (metal) melting. This behavior is present 
also in the modeling of Ghosh & McSween (1998), 



the scenarios (see Table |3|. 

In Table [3] we report the results obtained for dif- 
ferent values of Ai^ and post-sintering porosity, in 
particular the time at which the metal {rmet) melt- 
ing begins (T = 1213 K), the time of formation 
of core {T^ore), the time at which the silicate (t^j;) 
melting ends (T — 1850 K), the size of the surviv- 
ing chondritic crust, the size of the SL after 3 and 
5 Ma, the maximum degree of the silicate melting 
and the maximum temperature reached inside the 
asteroid after 5 Ma. The formation of the core (i.e. 
the separation of two melts) begins when the tem- 
perature reaches 1725 K in the region containing 



metals: 



indicates the time when all metals 



migrated toward the center. 

In the case of 1 vol.% porosity, the thickness of 
the chondritic crust ranges from 3 to 17 km, while 
the thickness of the solid layer ranges from 7 to 21 
km after 3 Ma and from 8 to 27 km after 5 Ma. In 
cases of porosity of 2 and 5 vol.%, the chondritic 
crust ranges from 4 to 18 km and 4 to 19 km, re- 
spectively. After 3 Ma the SL ranges from 8 to 21 
km and after 5 Ma the SL ranges from 9 to 30 km 
in both cases (see Table Isl). 

It is be noted that the temperature at the sur- 
face, while never reaching the solidus temperature 
of pure iron metal, can reach values as high as 1000 
K. It is also noteworthy that the SL thickness im- 
plies a ~5-10 km thick eucritic layer can already 
form between 3 and 5 Ma, in agreement with the 
dating of the oldest eucrites described by |Bizzarro 
eFaLl ((20051). 



Scenario 


Delay [Ma] 


Delay[r_Ai] 


NO 


0.00 


0.0 


Nl 


0.36 


0.5 


N2 


0.72 


1.0 


N3 


1.08 


1.5 


N4 


1.43 


2.0 


N5 


1.79 


2.5 


N6 


2.15 


3.0 



who however assumed this point to be fixed. In 
this work, instead, the location of this intersection 
moves in time depending on the characteristics of 



Table 2: Seven scenarios have been studied, where 
the delay-parameter /St^ is expressed in Ma and in 
half-lives of ^^k\. 
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5 Conclusions 

In this paper we investigated the thermal history 
of Vesta across the first 5 Ma of its history. We 
focused our analysis on the effects of time passed 
between the injection of ^^Al in the Solar Nebula 
and the formation of Vesta (which is assumed in- 
stantaneous), and consequently the content of ^^Al, 
and of the initial porosity on the subsequent evo- 
lution of the asteroid. Our model is based on the 
simultaneous study of the evolution of the internal 
structure (percolation of metals and separation of 
silicatic and metallic melts) and the thermal heat- 
ing of the asteroid, due to ^^Al and other short 
and long-lived radionuclides. We observed that the 
main source of energy is represented by ^^Al while 
the contribution of the other radionuclides is negli- 
gible. 

At the surface, the balance between thermal 
heating and the black-body radiation into space 
offers a more realistic picture of the thermal his- 
tory than a simple case in which the temperature 
of the surface is fixed to a constant value. We ex- 
plored several scenarios (N0-N6) characterized by 
different radiogenic strengths, expressed by a delay- 
parameter (Atd) ranging from to 2.16 Ma. We 
opted for an initial composition of about 75 vol.% 
of silicates and about 25 vol.% (similar to H or L 
chondrite compositions): the chosen values of the 
porosity are 1, 2 and 5 vol.%. In the scenarios in 
which the differentiation takes place a metallic core 
is always formed, its formation time ranging from 
about 0.5 to 3.5 Ma. We also observed the solidifi- 
cation of a surface layer whose maximum thickness 
is about 30 km. This solid layer is composed by an 
unmelted chondritic crust with its thickness reach- 
ing the maximum of about 20 km and the exact 
value strongly depending on the delay time (Ai^)- 
A ^5-10 km thick layer (which can be interpreted 
as the topmost part of the eucritic layer observed 
by Dawn) can already form beneath the chondritic 
crust between 3 and 5 Ma, in agreement with the 
dating of the oldest eucrites described by |Bizzarro] 
et al. (2005). The solidification of this layer is pos- 



sible even in absence of convection and volcanism, 
but only considering conduction of heat and radia- 
tion at the surface. 

It is interesting to note that the existence of 
a chondritic crust overlying a differentiated inte- 



based on analytic arguments, for the other two 
largest bodies in the asteroid belt, Ceres and Pal- 
las, but not for Vesta. Moreover, the independent 



results obtained by Tkalcec et al. (2013) using a 



2D finite difference code confirm our findings: while 
not explicitly discussed in their paper, their Figs. 4 
and S2 clearly show the existence of a surface layer 
of material that never experiences melting (i.e. our 
chondritic crust) whose thickness varies between 3 
and 7 km. A further comparison between the two 
models also reveals an agreement between the con- 
strains on the formation time of Vesta (G. Golabek, 
personal communication), i.e. Vesta's formation 
should have occurred preferably before 1 Ma and, 
more generally, no later than 1.5 Ma respect to the 
injection of ^^ Al in the Solar Nebula (which, in their 
work, coincides with CAIs crystallization). 

The survival of the chondritic crust and the rate 
of solidification of the underlying eucritic layer are 
particularly important to understand the geophys- 
ical and thermal evolution of Vesta as, across the 
temporal interval here investigated, the asteroid 
underwent a bombardment caused by the forma- 
tion of Jupiter (iTurrini et al.||2011l 120121). The Jo- 



vian Early Bombardment, in fact, causes a global 



erosion of the primordial crust of Vesta (Turrini 



submitted! and the chondritic crust, not currently 



observed by Dawn, plays an important role in pre- 
serving the eucritic layer. Moreover, the Jovian 
Early Bombardment can trigger local or large-scale 
effusive phenomena due to the excavation of craters 



and the formation of impact basins ( Turrini et al. 
2011 Turrini|submitted ) , thus affecting the cooling 



rior was speculated by Elkins-Tanton et al. (2011 ) 



history of Vesta. 

The complete melting of silicates is achieved if 
Vesta does not form later than 0.8 Ma after the 
injection of ^^Al in the Solar Nebula while it is 
achieved in a limited region of the mantle if the 
formation took place after 1 Ma: in the rest of the 
asteroid a 50 vol.% melting of silicates occurs. If 
the formation is completed no later than 1.4 Ma 
from ^®A1 injection the degree of melting is gener- 
ally lower than 50 vol.%. For delay times of more 
than 1.5 Ma the melting, when occurring, is limited 
in a small region of the asteroid. 

The formation of Vesta should not have occurred 
later than 1 Ma after ^^Al injection, if the crystal- 
lization of eucrites and diogenites are linked to a 
large degree (more than 50 vol.%) of silicate melt- 
ing ( Greenwood et al. 2005 ) . If eucrites and diogen- 
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ites can form from a partial melt ranging to 25-50 
vol.% the formation of Vesta could have occurred 
up to 1.4 Ma after ^^Al injection. 
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Tabic 3: Summary of scenarios. Wc report the time at which the metal melting occurs {Tmet, T = 1213 
K), the time of formation of the core (rcore), the time at which the silicate melting ends {tsu, T — 1850 
K), the chondritic crust thickness, the solid layer thickness after 3 and 5 Ma, the percentage of silicate 
melt and the maximum temperature reached after 5 Ma. 
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Figure 1: From left to right, from top to bottom: thermal history maps, for = 1 vol.%, of NO (a), Nl 
(b), N2 (c), N3 (d). Black isoline (1213 K) identifies the onset of metal melting; white (1425 K), black 
(1725 K) and magenta (1850 K) isolines identify the onset, the 50 vol.% and the complete melting of 
silicates, respectively. 
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Figure 2: From left to right, from top to bottom: thermal history maps, for = 2 vol.%, of NO (a), Nl 
(b), N2 (c), N3 (d). Black isoline (1213 K) identifies the onset of metal melting; white (1425 K), black 
(1725 K) and magenta (1850 K) isolines identify the onset, the 50 vol.% and the complete melting of 
silicates, respectively. 
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Figure 3: From left to right, from top to bottom: thermal history maps, for = 5 vol.%, of NO (a), Nl 
(b), N2 (c), N3 (d). Black isoline (1213 K) identifies the onset of metal melting; white (1425 K), black 
(1725 K) and magenta (1850 K) isolines identify the onset, the 50 vol.% and the complete melting of 
silicates, respectively. 
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Figure 4: From left to right, from top to bottom, plots are: temperature vs distance from center at 0.1 
Ma (a), at 0.5 Ma (b), at 1.5 Ma (c), at 3.0 Ma (d), at 5.0 Ma (e) and maximum temperature profile vs 
time (f), ior (f> — 1 vol.%. The horizontal lines represent: the windows for the melting of metals (green 
and red) and the of the silicates (cyan and magenta) . 
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Figure 5: From left to right, from top to bottom, plots are: temperature vs distance from center at 0.1 
Ma (a), at 0.5 Ma (b), at 1.5 Ma (c), at 3.0 Ma (d), at 5.0 Ma (e) and maximum temperature profile vs 
time (f), for (/) = 2 vol.%. The horizontal lines represent: the windows for the melting of metals (green 
and red) and the of the silicates (cyan and magenta) . 
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Figure 6: From left to right, from top to bottom, plots are: temperature vs distance from center at 0.1 
Ma (a), at 0.5 Ma (b), at 1.5 Ma (c), at 3.0 Ma (d), at 5.0 Ma (e) and maximum temperature profile vs 
time (f), for (jj — 5 vol.%. The horizontal lines represent: the windows for the melting of metals (green 
and red) and the of the silicates (cyan and magenta) . 
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